Controlling the transmission line shape of molecular t-stubs and potential 

thermoelectric applications 
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Asymmetric line shapes can occur in the transmission function describing electron transport in 
the vicinity of a minimum caused by quantum interference effects. Such asymmetry can be used to 
increase the thermoelectric efficiency of molecular junctions. So far, however, asymmetric line shapes 
have been only empirically found for just a few rather complex organic molecules where the origins of 
the line shapes relation to molecular structure were not resolved. In the present work we introduce a 
method to analyze the structure dependence of the asymmetry of interference dips from simple two 
site tight-binding models, where one site corresponds to a molecular n orbital of the wire and the 
other to an atomic p z orbital of a side group, which allows us to analytically characterize the peak 
shape in terms of just two parameters. We assess our scheme with first-principles electron transport 
calculations for a variety of t-stub molecules and also address their suitability for thermoelectric 
applications. 

PACS numbers: 



I. INTRODUCTION 

Despite increased interest in molecular electronics and 
tremendous advances in the theoretical description as 
well as experimental characterization of the conductance 
of single molecule junctions in the last two decades, there 
has been a shortage of novel ideas for molecular devices. 
Recently, two concepts emerged which received a lot of 
attention within the field: i) devices based on quantum 
interference (QI) effects 1 - - — and ii) theromelectric applica- 
tions (TA)£~— . QI devices are utilizing the wave nature 
of electrons and the resulting possibility for destructive 
interference of transmission amplitudes for the design of 
data storage elements^, transistors 2 - - — or logic gates 5 -, and 
QI is also the cause for the observed low conductance of a 
benzene contacted in the meta configuration as compared 
to the para and ortho configuration s 19 ' 20 . It has been 
rationalized by a variety of different physical pictures, 
such as phase shifts of transmission channels or interfer- 
ing spatial pathways^ - —. More recently, research in QI 
has been extended to aromatic molecules of increasing 
size 2 ^— as well as incoherent transport in the Coulomb 
blockade regim e 29 ! 30 and also led to proposals for the us- 
age of molecular interferometers in spintronic o 31 ' 32 . 

Thermoelectric materials on the other hand can con- 
vert thermal gradients to electric fields for power gener- 
ation or vice versa for cooling or heating which makes 
them useful as Peltier elements. The efficiency of a ther- 
moelectric material is measured by a dimensionlcss num- 
ber, the figure of merit ZT = S 2 TG/k, which includes 
the thermopower S 1 , the average temperature T, the elec- 
tronic conductance G, and the total thermal conductance 
k, which has contributions from electrons (K e ;) as well as 
from phonons (ttph)- hi Appendix A we show how these 
quantities are calculated. Materials with ZT ~ 1 are re- 
garded as good thermoelectrics, but ZT > 3 would be 



required to compete with conventional refrigerators or 
generators 3 - 3 -. The optimal thermoelectric material has a 
high thermopower and high electronic conductance but 
a low thermal conductance. One way to optimize ZT 
is thus to increase the thermopower. It is instructive to 
consider an approximate formula for the thermopower— 



T{nk B ) 2 d\n{T(E)) 



3e 



dE 



(1) 



where T{E) is the energy dependent electron transmis- 
sion function. This formula is valid at low temperatures 
and at energies not too close to transmission resonances. 
By inspection of Eqn. (1) it becomes obvious that a 
large thermopower can be achieved if the transmission 
function changes rapidly around the Fermi energy. This 
property therefore provides a link between QI and TA in 
the sense that QI effects can create steep declines in the 
transmission function for certain molecules. Although 
quite good thermo-electric characteristics have also been 
predicted for some systems with the QI induced mini- 
mum surrounded by two maxima symmetrically-^, ide- 
ally the minimum is closer in energy to one molecular 
resonance where the transmission has a maximum thus 
creating a very asymmetric transmission function, which 
naturally results in a steep transmission gradient in the 
respective energy range. In studies of mesoscopic quan- 
tum waveguides such asymmetric peak shapes caused by 
QI are usually referred to as Fano resonances ; 34 ' 35 which 
are not limited to electron transport but a quite general 
phenomenon in physics^ They occur due to the inter- 
ference between two wave amplitudes, one corresponding 
to a continuous state and the second to a localized state, 
and were originally observed in the line-shapes of inelas- 
tic electron scattering by helium atoms, where the con- 
tinuous state was defined by the background scattering 
process and the localized one by the process of autoion- 
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FIG. 1: (a) Transmission function, (b) thermopower, and (c) 
purely electronic contribution to the thermoelectric figure of 
merit ZT e i (i.e. the phonon contribution to the thermal con- 
ductance is neglected). The different curves are calculated 
from the general Fano-transmission function giving a sym- 
metric dip for q — (dashed red), an asymmetric Fano-shape 
for q = 1 (solid blue), and a symmetric peak for q — > oo 
(dash-dotted black). The temperature is T = 300 K and the 
broadening parameter is T = 0.2 eV. 



ization 3 -^. 

As an illustrative example of the effect of transmis- 
sion function asymmetry on S and ZT, we consider a 
model system with a transmission function described by 
the Fano function 



T(E) = 



(qT/2 + Ef 



(l + c2 2 )(£ 2 +r 2 /4)' 



Although the antiresonances found in mcsoscopic elec- 
tron waveguides with stub-tuners^— can be classified 
fg\ as Fano resonances, because there the direct transmis- 

sion supported by a continuum of states within a period- 
ically extended wire interferes with the states localized in 
the stub, the analogy does not hold for single molecule 
junctions where the continuum of the electrode states is 
interrupted by the junction which acts as a tunnelling 
barrier—. But even if the analogy is reduced to a sim- 
ilarity, asymmetric QI features (which at least resemble 
Fano resonances) have been observed also in the trans- 
mission functions describing electron transport through 
single molecule s 34 ' 44 . Moreover, the general principle 
that electron waves propagating from the left to the right 
electrode can interfere with localized states in the junc- 
tion can also apply in the coherent tunnelling regime. 
What has been lacking so far, however, is a system- 
ic) atic investigation of the dependence of the peak shape of 
QI induced transmission dips on characteristic features 
of the structure of the single molecule in the nanojunc- 
tion, which could match the level of understanding which 
has been achieved for mesoscopic waveguides. Recently, 
we established a relation between molecular structure 
and the occurrence or absence of QI transmission zeros, 
which could be formulated in terms of simple graphical 
rule s 45 ' 46 . In the current article we aim at developing 
a framework which allows to project the characteristic 
features of molecular wires with one side group, which 
are most similar to the mcsoscopic stub-tuner, onto sim- 
ple tight-binding (TB) models. Even for such relatively 
simple molecules there are structural elements such as 
the chemical nature of the side group or its exact po- 
sition on the wire contributing to the peak-shape in a 
way that is not self-evident. In our article we pursue 
a step-by-step approach to this problem where we first 
ask the question how the transmission peak shape and 
thermoelectric properties are related to the parameters 
in a two-site tight binding model in general mathemat- 
ical terms. Then we relate this model to the topology 
of the 7T states of conjugated organic molecules with side 
groups and where this process results in chemically stable 
molecular structures validate our conclusions with first 
(2) principle electron transport calculations. 



The defining feature of a Fano resonance is that the 
resonance profile of the scattering cross section can be 
brought into this form. In Fig. [T] (a) we plot the trans- 
mission function T(E) in Eqn. © for three different val- 
ues of the shape parameter q giving rise to a symmetric 
dip (q = 0), an asymmetric Fano-shape (q = 1) and a 
symmetric transmission peak (q — >• oo). Panels (b) and 
(c) in Fig. [T] show the corresponding thermopower and 
ZT, respectively. Clearly the asymmetric transmission 
function gives the largest thermopower and ZT. At lower 
temperature and smaller T-values, the transmission peak 
will result in significantly higher S and ZT, but the asym- 
metric Fano-shape will always give larger values than the 
symmetric transmission dip. 



II. TWO-SITE TIGHT BINDING MODEL 

The TB-modcl, which we introduce in the following 
aims at a separation of the structural aspects governing 
the peak shape of the interference dips, where analytical 
expressions can be derived for the transmission function 
in dependence on just two parameters. For achieving 
this aim we make a few simplifying assumptions. First, 
we assume that only the ir electrons are relevant in the 
energy range of interest, which allows us to describe the 
molecules within a TB model with only one p z atomic or- 
bital (AO) for each carbon- or hetero-atom. In a second 
step, we look at molecular orbitals (MOs) of the t-stub 
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FIG. 2: Analytical TB results for a two-site TB model (f) 
according to Eqns. (l)-(3) with 7 = 0.5 eV and £0 = 0.0 eV. 
(a) Transmission function T as a function of the electron in- 
cident energy E (with ei = 0.5 eV) and (b) AE as a function 
of the side group AO energy ei with the MO-AO coupling 
|*| = 0.1 eV (red), 0.5 eV (black) and 1.0 eV (green), (c) 
T(E) (with \t\ = 0.5 eV) and (d) AE(\t\) for ei = 0.0 eV 
(red), 0.5 eV (black) and 1.5 eV (green). Panel (e) shows the 
purely electronic thermoelectric figure of merit ZT e i calcu- 
lated from the transmission functions in panel (c). 



wires without the side group attached, which can be ob- 
tained by diagonalizing the TB-Hamiltonian for a chain 
of carbon atoms with a certain length and from them pick 
the single MO which is closest in energy to the p z AO 
or fragment orbital of the side group. If we now assume 
a two-site single-particle picture where the transmission 
function for electrons incident on a molecular junction 
with an energy E is dominated by this MO and its inter- 
action with the side group AO, it can be written as 



T{E) 



4 7 2 



(E-e 



E-e, 



-4 7 2 



(3) 



as derived in Appendix B, with 7 being the lead cou- 
pling within a wide band approximation, Eq the wire MO 
energy, E\ the side group AO energy and t the coupling 
between them. For an illustration see Fig. [5] (f). As ex- 
plained above the thermoelectric properties of a molecu- 



lar junction can be optimized if the transmission function 
rapidly changes from a minimum to a maximum. As a 
simple measure of how rapidly the transmission function 
changes we consider the energetic difference AE between 
the transmission minimum and the closest maximum. It 
follows from Eq. ([3|) that the transmission is zero at en- 
ergy Eq = £ir^ while for transmission peaks with T = 1 



the quadratic equation E 



E — Ae 



has to be solved 



where we made the substitution Ae — E\ — Eq and set 
£0 = without loss of generality for our following argu- 
ments. The solutions to this equation E\^ and the norm 
of their differences with the energy of the transmission 
zero Eq can then be obtained as 

Ae 
2 

and 

Ae 
2 



■Ei, 2 



± 1 -^/Ae T 



4i 2 



I -Ei, 2 — E \ — I — 



4i 2 | 



(4) 



(5) 



where it can be seen that only for finite Ae there will al- 
ways be one solution for \Ei^ — E \ which will be closer to 
zero than the other. In other words only under this con- 
dition will the energy difference of the two transmission 
maxima from the one transmission minimum be asym- 
metric for the two-site model, while for Ae = we get 
|Ei,2 — Eq\ = ±£ and the peak shape around the minimum 
in T{E) will therefore be symmetric. In the following we 
will refer to the smaller of the two solutions for \E\.2~Eq\ 
as AE, which we use as a measure indicative of thermo- 
electric efficiency where we expect the thcrmopower to 
be high for small AE. We also find in practice that in 
most cases a small value for AE indicates a high level 
of asymmetry in the peak shape around the minimum of 
the transmission function but this argument cannot be 
made in general, where a ratio of the two solutions should 
be used instead. Another conclusion at this point is that 
|^i,2— Eo| is completely independent of the lead-coupling 
7. This independence of rapid transmission changes from 
lead-coupling broadenings is a general advantage of uti- 
lizing QI since both transmission peaks and transmission 
minima are defined entirely by the molecular topology. 
Temperature effects can, however, smear out very sharp 
features when ksT > AE. 

In Fig. [2] we illustrate explicitly for the two-site model 
how the shape of the transmission function changes with 
the size of t (panel (a)) and Ae (panel (c)), where the 
peak shapes are found to be the more asymmetric, the 
larger Ae is and the smaller t. The variation of AE with 
Ae and t are shown in panels (b) and (d). We note that 
these two parameters correspond to the length and height 
of the tunnelling barrier introduced in Rcfsi 40 ' 41 for a sep- 
aration of the sidearm in mesoscopic waveguide t-stubs, 
where the authors also found that distinctly asymmetric 
line shapes only start to appear for weak coupling and 
large barrier heights. 

Figure [5] (e) shows the purely electronic ZT e i (phonon 
contribution to thermal conductance is neglected) cor- 
responding to the transmission functions in panel (c). 
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Clearly the maximum ZT e i values increase when the 
transmission function becomes more asymmetric, in cor- 
respondence with the findings in Fig. [TJ While very large 
ZT e i values have been predicted for symmetric transmis- 
sion functions 1 ^, the results in Fig. [2](e) show that asym- 
metric transmission functions are even more promising 
for a thermoelectric purpose. 



III. MOLECULAR ORBITAL ANALYSIS 

In order for our two-site model to be useful for the 
analysis of t-stub molecular wires we need to be able to 
deduce the parameters t and As from the full topology 
of the wires p z AOs including the side group, which docs 
not always have to be another carbon-site. For this we 
start from an AO-TB-model where the onsite energy of 
each p z orbital of the carbon chain eq is set to zero and 
the couplings between neighbouring sites a arc defined 
as -2.9 eV. From a diagonalization of the chain without 
the side group, where we refer to Appendix C for the full 
mathematical details, we obtain its MO energies and the 
orbital weight on each chain site for each MO. Within 
this scheme we can vary the side group AO energy e sg 
and an initial value for the MO-AO coupling j3 freely 
so that it reflects the side groups chemical nature. This 
initial value for (3 is then multiplied by the orbital weight 
at the chain site to which the side group is bonded in 
order to derive the coupling parameter t sg . The relevant 
MO is typically the one closest in energy to e sg with a 
finite value for t sg . The two site model finally contains 
this latter MO (with energy Emo) an d the side group AO, 
where the lead coupling 7 in the two-site model is derived 
from the wide band lead coupling T used in the full AO- 
TB model multiplied by the squared orbital weight on the 
terminal sites in the chain. In Figs.[3]and [4] we illustrate 
how this works in practice and which kind of information 
can be obtained from such an analysis. At this stage we 
distinguish the chemical nature of the side groups only 
by their onsite energies and define (3 — — 2.4eV for all 
the calculations illustrated in the two figures. But even 
with (3 fixed to a given value, t sg does vary with the MO 
topologies as discussed below. The side groups in our 
AO-TB model are always single AOs which for multi- 
atom groups such as CH2 and NO2 represent fragment 
orbitals (FOs) where we explain in Appendix D how to 
obtain from DFT calculations the onsite energy of the 
FO most relevant to the QI effect. 

We want to stress that the AO-TB model we use in 
this section is only meant to investigate topological ef- 
fects in a meanfield description. As we point out in the 
following, some of the topologies we study here would 
represent radicals or diradicals if they are translated into 
chemical structures and for a realistic description of their 
behaviour in a single molecule transport junction a multi- 
determinant approach would be neede d 47 i 48 . This is, 
however, not our focus in this article where the AO-TB 
model is meant to provide a bridge between the 2-site 



analysis in the previous section and the NEGF-DFT cal- 
culations in the next section. For the latter we limit our 
discussion only to junctions which are chemically stable 
and well described within a single-determinant approach. 

Fig. [3] focuses on molecular wires containing a side 
group X with s sg = — 0.5eV simulating a CH2 group. 
This is attached to a chain with nine atoms (in the fol- 
lowing referred to as C9) at the fifth (C9-5, Fig. (3k) and 
fourth (C9-4, Fig. [3p) site, respectively, and to a chain 
with eight atoms (in the following referred to as C8) at 
the fourth site ((C8-4, Fig. |3b). In the insets of Fig. 
we show the molecular topologies, where only C9-5 rep- 
resents a chemically stable system, while C8-4 would be 
a radical and C9-4 even a diradical if the side group is 
linked to the chain by a double bond. Nevertheless it is 
instructive to compare the three scenarios on a TB level 
in order to show that while there is no topological limit 
preventing transmission peak asymmetry for molecular 
chains with CH2 side groups, in practice such molecules 
will always have rather symmetric transmission peaks 
due to a selection criterion defined by chemical stabil- 
ity. 

The solid black lines in the transmission functions in 
the upper panels of Fig. [3] are the results of full TB AO 
calculations including the side group, while the dashed 
red and blue lines come from 2-site TB MO models con- 
sidering only the side group AO and the HOMO (red) 
or HOMO-1 (blue) of the unsubstituted chain, respec- 
tively. For C9-5 the QI feature in T(E) is quite sym- 
metric and the black line well reproduced by the red line 
meaning that just considering the HOMO and the side 
group AO is sufficient for describing the line shape. The 
green line in Fig.[3£i shows the TB AO result without the 
side group, where it can be clearly seen that the peak in 
T{E) corresponding to the HOMO splits into two when it 
hybridizes with the side group AO, whereas the other two 
visible peaks corresponding to the HOMO-1 and LUMO 
remain unaffected. The situation is quite different for 
C9-4 (Fig. [5Jd), which exhibits a distinctly asymmetric 
QI feature and the red line does not reproduce the black 
one at all and instead represents a single site Lorentzian 
peak at the HOMO energy. The blue line on the other 
hand shows very good correspondence with the TB AO 
transmission, which means that the side group AO inter- 
acts mainly with the HOMO-1 and not with the HOMO. 

Our findings can be understood in terms of the topolo- 
gies of the MOs plotted in the lower panels of Fig. [3l 
where the HOMO has a high weight on site 5 but none 
on site 4 and in contrast the HOMO-1 has its largest 
weight on site 4 and none on site 5. If Emo defines the 
energy of the MO closest to £ sg with a relatively high 
coupling parameter t sg , these differences in MO topolo- 
gies also mean that for C9-5 |Ae| = \e S9 — £mo\ = 0.5 eV 
is substantially smaller than the value for C9-4 (1.49 eV), 
which explains why in the latter case the QI line shape is 
much more asymmetric according to the general trends 
shown in Fig. [5J We note that due to the symmetry prop- 
erties of the MOs relevant to our analysis C9-3 gives a 
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FIG. 3: Transmission functions (upper panels) and MO-AO hybridisation schemes (lower panels) for (a) C9-5, (b) C9-4 and 
(c) C8-4 type molecular stubs with a side group AO X (e S9 =-0.5 eV), where the respective AO topologies are given as insets 
in the transmission figures. The labels 5 and 4 correspond to the atom number in the carbon chain to which the side group 
is coupled, where we indicate the numbering we use for the atoms in the carbon chain in the lower part of panel (a). The 
line types for T(E) correspond to full TB AO calculations with (black) and without (green) the side group, and a 2-site MO 
model with the side group AO and the HOMO (ehomo=0 eV, dashed red) or the HOMO-1 (effOAfO-i=-L8 eV, dashed blue), 
respectively. The Fermi level E_f is assumed to be identical to ehomo of C9 (without any side groups). The lower panels also 
show the topologies of the HOMO and HOMO-1 where the respective position of the side group is indicated by red empty 
spheres, where depending on |Ae| and \t sg \ different hybridization patterns are observed. 




FIG. 4: Transmission functions (upper panels) and MO-AO hybridisation schemes (lower panels) for the AO topologies (which 
are shown as insets in the transmission figures) (a) C9-5 with a side group AO Y (£ sg =-2.5 eV), (b) C9-5 and (c) C8-4 with a 
side group AO Z (e S9 =1.5 eV). The low on-site energy of AO Y is energetically closest to the HOMO-1 (a), but since this does 
not have any orbital weight on the central carbon atom 5 it does not interact with the side group and produces a Lorentzian 
transmission function (dashed blue). A reasonable description of the transmission minimum is obtained from the HOMO-2 
and the side group (dashed green). In a similar way, AO Z has a side group energy close to the C9 LUMO (b), but since it 
has no orbital weight on chain site 5 it does not couple to the side group (dashed blue), while the interaction with the HOMO 
(dashed red) becomes the dominant one with a good qualitative description of the full transmission (solid black). For C8-4-Z 
(c) the side group is again close to the LUMO, which does have a finite orbital weight on chain site 4 and a 2-site model with 
the LUMO and the side group (dashed blue) closely reproduces the full TB-AO result (solid black). 
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transmission function similar to that of C9-5, and C9-2 is 
similar to C9-4, which we checked but do not show here. 
If the carbon chain has one site less (C8, Fig. |3t) no such 
dependency on the position of the side group can be iden- 
tified simply because the HOMO has a finite localisation 
on all sites 3-6, and therefore |Ae| = 0.5 eV for all cases. 
Due to the spread out HOMO localisation pattern, how- 
ever, t sg = 0.73 eV for C8-4, which explains why it shows 
a somewhat more asymmetric line shape than C9-5 in 
Fig. EJi, where t = -1.07 eV. The topologies in Fig. [3] 
cover all the inequivalent possibilities for wires with nine 
or eight carbon atoms in the chain with a single CH2 
side group. We found that only one of them is chemi- 
cally stable (C9-5) and this one has a symmetric trans- 
mission minimum. Since there is no reason to assume 
that shorter or longer wires would result in MO topolo- 
gies markedly different from those investigated, one has 
to conclude that such t-stub molecules with a CH2 side 
group are unlikely to exhibit asymmetric transmission 
shapes in general. 

Fig. g] now addresses the question what happens if we 
choose different side group AOs Y and Z with onsite en- 
ergies e sg — — 2.5 and 1.5 eV, respectively for simulating 
the effect of O and NO 2, while we keep the same value 
for p as in the study on the CH2 substituted wires with 
side group orbital X above. Interestingly, although the 
position of the minimum in T{E) is different (since it is 
always identical to e sg for molecular wires with one side 
group^i) , it can be already said by visual inspection that 
there is a close correspondence of the peak shape patterns 
encountered in Fig. [3] and those found in Fig. g] Figs. 0Jt 
and b, which mimic C9-5-0 and C9-5-N02 wires, respec- 
tively, are very similar to the one mimicking the diradical 
C9-4-CH2 in Fig. [3b>. This is because they exhibit again 
quite high values for |Ae|, namely 0.91 (with respect to 
the HOMO-2) for C9-5-Y and 1.5 with respect to the 
HOMO for C9-5-Z. C8-4-Z (Fig. gb) on the other hand 
resembles C8-4-X (Fig. [3b). In both cases |Ae| = 0.5 and 
t sg = 0.73 eV, although the relevant MO is the HOMO 
in Fig. [3b and the LUMO in Fig. gb. 



IV. FIRST-PRINCIPLES ELECTRON 
TRANSPORT CALCULATIONS 

As a next step we test our TB predictions and their 
interpretation on the basis of 2-site MO models by per- 
forming DFT based electron transport calculations^ for 
the transmission function of the four molecules shown in 
the insets of Fig. [5] As discussed in the previous section 
some topologies cannot be realised chemically, since they 
would be chemically unstable. The molecules C9-CH2- 
5, C9-0-5 and C8-N0 2 -4 in Fig. [5] on the other hand 
are legitimate targets for organic synthesis. While this 
is not true for the radical C9-N02-5, there is still some 
justification for studying the respective junction with a 
DFT approach, since the nitro group can stabilize radi- 
cals at least as transition states in conjugated ir systems 




FIG. 5: Transmission functions calculated for the molecular 
junctions shown in the inset from full NEGF-DFT calcula- 
tions^ (black lines), and from 2-site MO models (red lines) 
with the parameters |Ae| and \t sg \ as listed in Table [I] and 
obtained directly from DFT. For details see text. 



via the mesomeric effect. We also derive the TB param- 
eters for projecting the transport characteristics of these 
molecules onto 2-site models from the same calculations. 
Our procedure for the latter has the following steps: i) 
The local eigenstates on the side group orbital is found 
from a subdiagonalization of part of the full Hamiltonian 
corresponding to the LCAO basis orbitals on the side 
group atom. This yields information of the side group 
energy and the coupling to the nearest carbon p z orbital 
in the chain (see Appendices C and D for more details). 

ii) For projecting the topology of the molecular tt states 
without the side group we use the same scheme as with 
the AO-TB-models above, where the lead coupling 7 now 
is used as a fitting parameter in order to ensure that the 
broadening in the transmission functions calculated from 
the 2 site-MO model resembles the DFT transmission. 

iii) Finally we allow for a rigid energy shift of the 2 site- 
MO transmission function in order to account for Fermi 
level alignment in the DFT calculations^— . We show 
T(E) for all four molecules as calculated from NEGF- 
DFT (black lines) and from the corresponding 2 site-MO 
model calculations (red lines) in Fig. [5] and list the ob- 
tained model parameters |Ae| and \t sg \ as well as the 
values for Eq and \AE\ obtained from the black curves 
in TableU 

The peak shapes in Fig. [5]which are now based on DFT 
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molecule 


C9-CH 2 -5 


C9-0-4 


C8-NO2-4 


C9-N0 2 -5 


|A e | 


0.40 


0.91 


0.61 


1.60 




1.03 


1.12 


0.32 


0.40 


E 


0.26 


-1.44 


1.61 


1.72 


\AE\ 


0.89 


0.65 


0.40 


0.20 


S(E) 


-170 (0.21) 


184 (-1.39) 


199 (1.65) 


349 (1.74) 


ZT el (E) 


0.82 (0.19) 


0.93 (-1.37) 


1.00 (1.67) 


3.00 (1.77) 


ZT (E) 


0.003 (0.19) 


0.01 (-1.37) 


0.03 (1.67) 


0.48 (1.77) 



TABLE I: 2-site MO model parameters |Ae| and \t sg \ used 
for the red lines in Fig. [5] The energy of the transmission 
zero Eo and its difference to that of the the nearest transmis- 
sion maximum | AE\ are evaluated from the DFT transmission 
curves (black lines in Fig. 0. The maximum values for the 
thermopower S and the electronic contribution to the figure 
of merit ZT e i(E) (obtained at the corresponding energies in 
parantheses) and the full figure of merit ZT (including phonon 
contributions and calculated at the same energy) are listed 
in the bottom rows. All energies and the coupling strength 
\t sg \ are given in eV, while S is defined in fiV/K, and ZT e i 
and ZT are dimensionless and calculated at the temperature 
T = 300 K. 



calculations are very similar to those obtained from the 
TB models in Figs. [3] and 2J Now, however, the hopping 
element f3 has realistic values of -2.4, -2.5 and -1.5 eV 
for CH2, O and NO2, respectively, which naturally has 
an effect on \t sg \ as listed in Table HI The smallest value 
for \AE\ is found for C9-N0 2 -5 with the largest |Aer| 
and a rather small \t sg \, while C9-CHV5 has the largest 
|A2£| with the smallest |Ae| and a large \t sg \. In both 
cases the side group couples mostly to the HOMO of the 
unsubstituted molecule. But because of the difference 
in the onsite energies of the side group, |Ae| differs by 
1.2 eV between the two, while the weaker coupling of 
NO2 to the conjugated ir electrons of the wire is reflected 
in a reduction of a factor of two in |t Sff |. Both effects 
combine to explain the factor of four in \AE\ when these 
two systems are compared. C9-0-5 and C8-N0 2 -4, where 
the side group couples strongest to the HOMO-2 and the 
LUMO, respectively, are in between these two extremes 
in terms of |A£J|, where the low coupling parameter \t sg \ 
for the nitro group seems to dominate over the larger 
value of |Ae| found for the oxygen side group. 

Finally we want to address the suitability of the four 
molecules in Fig. [5] for thermoelectric devices, where we 
calculate the thermopower S and the figure of merit ZT 
explicitly from DFT (for a description of the theoretical 
framework we use for that we refer to Appendix A and 
Refspii~— ). We calculate both the purely electronic ZT e i, 
which is completely determined by the electron transmis- 
sion function assuming a vanishing phonon thermal con- 
ductance, K p h, and the full ZT with a constant value of 
Kph = 50pW/K as it was determined experimentally for 
quite similar alkanes with thiol anchors^. 

Tabic [T] lists the calculated thermoelectric properties 
for the four molecules in Fig. [5] Notably the qualitative 



ranking in S and ZT reversely corresponds to that of 
|A£J|, reflecting that the more asymmetric T(E) curves 
result in higher values for thermopower and figure of 
merit. This remains true and is at the highest end even 
emphasized when the contributions of phonons are con- 
sidered. The two drawbacks for C9-N02-5 in terms of 
applications are the rather high value of Eq and its in- 
stability as a radical. These results demonstrate that for 
an efficient and robust way of implementing QI induced 
transmission minima for thermoelectric applications an 
asymmetric transmission function with a transmission 
peak located close to the minimum in a chemically stable 
structure is needed. Since our analysis covers a large part 
of the chemical possibilities for t-stub type molecules, i.e. 
simple wires with one side group, one has to conclude 
that more complex probably aromatic molecules need to 
be investigated systematically with regard to their suit- 
ability for thermoelectric devices. 



V. SUMMARY AND CONCLUSIONS 

In summary, we presented a scheme for the systematic 
investigation of the dependence of the symmetry of the 
transmission function around the energy of quantum in- 
terference induced minima on the molecular structure in 
single molecule junctions. For this purpose we mapped 
the molecular orbitals of molecular wires with side groups 
onto simple two-site tight binding models where it was 
found that the shape of the transmission function just de- 
pends on two parameters which provides a strong link to 
previous findings for t-stub tuners in mesoscopic waveg- 
uides. We applied the mapping procedure on a series of 
molecular wires which vary in their number of carbon 
atoms as well as the chemical nature and bonding posi- 
tion of side groups. We predicted from our scheme and 
verified with density functional theory calculations that 
the peak shape is mostly determined by the coupling be- 
tween the side group and the molecular chain (where a 
small coupling results in high asymmetry) and the en- 
ergy difference of the side group atomic orbital and the 
molecular orbital closest to it (where a large difference 
results in high asymmetry). The role of the topology of 
this molecular orbital in producing quite different peak 
shapes for very similar molecules could be explained from 
our model and was further verified by DFT transport 
calculations. Asymmetric peak shapes are highly desir- 
able in terms of device applications, because they promise 
high thermoelectric efficiency, and the projection scheme 
we presented for the systematic study of peak shapes 
therefore provides a route towards the chemical engineer- 
ing of thermoelectric devices. 
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Appendix A: Figure of merit - ZT 

The efficiency of a thermoelectric material can be char- 
acterized by the dimcnsionless figure of merit, ZT, given 
by 



ZT 



S 2 GT 



(Al) 



where S is the Seebeck coefficient, G the electronic con- 
ductance, T the temperature, and n p h and n e i are the 
phonon- and electron contributions to the thermal con- 
ductance, respectively. Except the phonon contribution 
to the thermal conductance n, p h, all quantities can be 
obtained from the electron transmission function. 
The linear electronic conductance is 



where we have introduced the function L m (ii): 

d.f(E,fi) 



2 r°° 

Uri = 7 &ET{E){E-nY 
h J-oo 



dE 



(A3) 

The thermopower, S, is defined as the voltage required 
to counterbalance the electronic current flowing due to a 
(small) temperature gradient across the junction 



S 



lim — — 

AT^O AT 



1=0 



(A4) 



while the purely electronic contribution to the ther- 
mopower can be written compactly as^ 4 - 



S(fi,T) 



eT L (/x) ' 



(A5) 



There is in principle also a so-called phonon drag con- 
tribution to the thermopower due to electron-phonon in- 
teractions. This term is usually small and we shall not 
consider it here. 

Finally, the electronic contribution to the thermal con- 
ductance is given by^ 4 - 



u ' T \ L (n) 



It is useful to rewrite ZT as 
S 2 GT 



ZT 



K,:l 



K e l 



ZT ei 



K e l 



Iph + K el 



(A6) 



(A7) 



where ZT e i is the purely electronic part which equals the 
total value for ZT, when n e i 3> K. p h- In recent theo- 
retical articles on thermoelectrics in molecular junctions 
very large values for ZT e i have been reported for molecu- 
lar junctions with quantum interference effects and tran- 
mission node s 12 ' 13 . However, it is not really justified to 
neglect n p h, especially not when the electronic transmis- 
sion is very small, as is the case around transmission ze- 
ros, because in this case K e i will also be very small, and 
the contribution from phonons may be very important. 

While calculations of the phonon thermal conductance 
in principle can be performed at a DFT level of the- 
ory, the calculations are quite demanding and there are 
only few reports on first-principles calculations of K p h for 
molecular junction s 14 ' 18 . The main focus in our present 
work is on electronic properties and the ZT values re- 
ported in the main text are thus calculated from Eq. (|A7[) 
assuming a constant value of K p h = 50 pW/K as deter- 
mined experimentally for thiol-bonded alkanes^ 3 -, which 
are quite similar to the molecules considered in our arti- 
cle. 



Appendix B: Derivation of the transmission function 
for two TB orbitals in the wide band limit 



Assuming a single-particle picture, the transmission 
function for electrons incident on a molecular junction 
with an energy E can be defined as 



T(E) = Tr[Gr L G t r fi ](£) 



(Bl) 



where G = (EI — H mo i — — S^) -1 is the Green's 
function matrix of the contacted molecule, I is the iden- 
tity matrix, T,l/r is the self-energy due to the left/right 



lead, and T 



L/R 



L/R 



In our AO-TB mod- 



els the Hamiltonian describing the molecule is given in 
terms of a basis consisting of localized atomic-like or- 
bitals, <f>\, . . . ,4>m , j where only two are coupled to the 
leads, which allows us to write the transmission function 



as 



T(E) = 4y(E) 2 \G 1N (E)f 



(B2) 



The matrix element Gin(E) can be obtained using 
Cramer's rule 



G 1N (E) 



C 1N (EI-H 



det(E - H mol - S L - 



(B3) 



where C\n{EI — H mo i) is the (IN) co-factor of (EI — 
H mo i) defined as the determinant of the matrix ob- 
tained by removing the 1st row and iVth column from 
(EI - H mo i - Si - Er)) and multiplying it by (~1) 1+N . 
Since we assume that only orbitals <f>\ and 4>m couple to 
the leads, the removal of the 1st row and Ath column 
completely removes £z,,i? in the co-factor. We imply a 
wide band limit, where the lead coupling strength, 7, 
becomes energy independent and S 



t 

L/R 



-vy. 
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For the two-site MO models in the main text there is 
only one site with an onsite energy e connected to both 
leads (N — 1) and coupled to a side group (ei) with a 
hopping parameter t. For such a system we find 



which results in 



Cn(E — H mo i) — E — ei 
det(£ - H mol -H L - Sfl) 



and 



E - e - 2-i7 

-i £ - ei 



(B4) 
(^5) 



T(E) 



(E - £l ) 2 4 7 2 



4 7 2 



((i?-eo)(i?- e i)-t 2 ) 2 + 4 7 2 (i ; ;-e 1 ) 2 

I 



(B6) 



Appendix C: Derivation of the MO projection 
scheme 



In order to address the explicit dependence of the side 
group on-site energies and coupling parameters we first 
write the full Hamiltonian of the molecule (in a basis of 
p z orbitals) as 



H — H o + H sq , 



(CI) 



where H describes the linear carbon chain without the 
side group having the structure 



( e c a 
a e c a 



H 



\ o o 







(C2) 



0/ 



Here e c is the carbon on-site energy and a is the nearest 
neighbour hopping parameter. The last row and column 
correspond to the single basis function on the side group. 
The side group Hamiltonian, H sg has the form 



/ 



H, 



\0 



a o 



o \ 

p 





(C3) 



; »g 



We now proceed by diagonalizing H to find the molecu- 
lar orbitals (MOs) of the molecule without the side group. 
The i'th MO has the vector form 



i^ ) ) = [c 1 l) , 4°, 



(C4) 



(i) 

assuming there are N atoms in the chain, and where Cj 
is the coefficient corresponding to the orbital weight on 
site j. The last '0' is at the side group orbital where the 
MO by definition has no weight. We now make the ap- 
proximation, to consider only a limited number of MOs 
and neglect all the others. In the simplest case, we con- 
sider only the HOMO, \ip H ), and proceed to construct 
a model Hamiltonian including only the HOMO and the 
side group orbital. In this case our basis consist of two 
orbitals: 



l<M = \^ H ) 

\<t> 2 )=e\r 9 ) 



[c? 
[0, 0, 



i L 2 ' 



...,4,0] (C5) 
,0, 1], (C6) 



where the first one is a MO of the molecule without side 
group and the last one is the side group orbital. The 
matrix elements of the model Hamiltonian are given by 



hi 



biWj) = (^\H \<p 3 ) 



\H M 



(C7) 



Since (0i|-ff o |<M = e H (the energy of the HOMO), 
(0 2 | W2) = 0, (0i|iT ifl |0i) = 0, and {<t> 2 \H sg \<t> 2 ) = e sg , 
we arrive at 



h = 



£H tsg 
tsg £sg 



where the coupling parameter t sg is given by 



h\H + H s 



(C8) 



(C9) 



and Cj is the MO orbital weight on the site, which cou- 
ples to the side group. We thus observe that the effective 
coupling, t sg , to the side group in the model Hamiltonian 
depends on the real coupling parameter (3 and the orbital 
weight. This described procedure where only one MO is 
considered is applied for all the molecules in the main 
text. The transmission through the model system is cal- 
culated in the wide band limit with constant self-energies 
assuming that only one MO couples to both electrodes. 
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The elcctrodc-sclf energy matrices thus have the form 

and the transmission function is calculated using the 
standard non-equilibrium Green's function (NEGF) 
methods also employed in conjunction with density func- 
tional theory (DFT) in the main text. 

The value of 7 is obtained from the lead coupling T 
used in the full AO-TB model as 

il = r L | c f | 2 (en) 

1R = T R \4\ 2 . (C12) 

The reason why the orbital weight is squared is that the 
lead self-energy, is calculated from a surface Green's func- 
tion, g as T,l.r = v£ R gVL,n, where the effective cou- 
pling, V, between the leads and the molecule can be cal- 
culated for the single HOMO orbital in the same way as 
its effective coupling to the side group shown above. In 
all our calculations there are no differences between left 
and right coupling since all the considered orbitals have 
the same squared weight on both terminal sites of the 
chain. 



Appendix D: Computation of on-site energies 

We calculate the side group on-site energy from the full 
Hamiltonian matrix, H, describing the molecule and the 



Au electrodes. We project onto the subspace spanned by 
the basis functions of the side groups 



h — P HP 

L1 sg A sg ±±x - sg> 



where P sg has diagonal elements on the indices of the side 
group basis functions, and zeros elsewhere and similarly 
obtain also a side group overlap matrix, s sg . We then 
diagonalize s^hgg to find the side group energies and 
eigenstates, where for each side group eigenstate \ipi), 
we calculate a coupling matrix element to the closest p z 
orbital (assumed to be on atom n) in the main chain as 



= (p n z \mi), (di) 



with \p") a vector which is one on the index correspond- 
ing to the p z basis orbital on atom n and zero in all 
other coordinates. Only the side group eigenstates with 
7r symmetry will result in a large coupling tf. We there- 
fore include only the tt state which is closest to the Fermi 
level in our 2 site MO model and take its energy for the 
parameter e sg . 
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